Reference: Biol. Bull. 204: 10-20. (February 2003) 
© 2003 Marine Biological Laboratory 


Molecular Genetic Evidence that Dinoflagellates 
Belonging to the Genus Symbiodinium 
Freudenthal Are Haploid 


SCOTT R. SANTOS* AND MARY ALICE COFFROTH? 

DepurUueiit of Biological Science, State University of New }ork at Buffalo, 
Buffalo, New Yoi'k 14260-1300 


Abstract. Microscopic and cytological evidence suggest 
that many dinotlagellates possess a haploid nuclear phase. 
However, the ploidy of a number of dinoflagellates remains 
unknown, and molecular genetic support for haploidy in this 
group has been lacking. To elucidate the ploidy of symbiotic 
dinoflagellates belonging to the genus Syjuhiocliniuin, we 
used five polymorphic microsatellites to examine popular 
tions harbored by the Caribbean gorgonians Plexaiira kuna 
and Pseudopterogorgia elisahetliae: we also studied a series 
of Sytnbioduiiinu cultures. In 690 out of 728 Synibiodinium 
samples in liospite (95% ol the cases) and in all 45 Syni¬ 
biodinium cultures, only a single allele was recovered per 
locus. Statistical testing of the Symbiodinium populations 
harbored by P. elisabetliae revealed that the observed ge¬ 
notype frequencies deviate significantly from those ex¬ 
pected under Hardy-Weinberg equilibrium. Taken together, 
our results confirm that, in the vegetative life stage, mem¬ 
bers of Symbiodinium, both cultured and in liospite, are 
haploid. Furthermore, based on the phylogenetics of the 
dinoflagellates, haploidy in vegetative cells appears to be an 
ancestral trait that extends to all 2000 extant species of these 
important unicellular protists. 

Introduction 

The ploidy of an organism can significantly affect ge¬ 
nome evolution. For example, diploids caiTy twice as much 
DNA as haploids and may be expected to accumulate new 
beneficial mutations at a higher rate (Paquin and Adams, 
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1983). The genome of diploids may also evolve rapidly 
because they carry more than a single copy of an allele. 
These “extra" alleles, over time, may evolve new functions 
while “old" alleles continue to perform their original func¬ 
tions (Lewis and Wolpert, 1979). Haploids, on the other 
hand, tend to have deleterious mutations purged more rap¬ 
idly from the population since they are not masked (Hughes 
and Otto, 1999). Furthermore, knowledge of ploidy is es¬ 
sential to the inteiprctation and understanding of population 
genetic data. Given the importance of ploidy to genome 
evolution and population genetics, it is surprising that ques¬ 
tions pertaining to it remain for a number of organisms. 

Dinotlagellates are a diverse and ecologically important 
group of unicellular protists. For example, some species arc 
major photosynthetic or heterotrophic components of the 
plankton, and others are considered to be the causative 
agents of fish kills. Microscopic and cytological evidence 
from the species examined to date suggests that dinollagel- 
lates (with the exception of Noctiliica .spp.) possess a veg¬ 
etative haploid nuclear phase (Phester and Anderson, 1987; 
Coats, 2002). However, ploidy has not been explicitly de¬ 
termined for a number of dinotlagellates, including the 
important genus Symbiodinium Freudenthal (Taylor, 1974). 
Members of Synibiodinium, commonly referred to as zoo- 
xanthellae, arc intra- or intercellular symbionts of manne 
invertebrates, including foraminiferans, sponges, cnidar- 
ians, and molluscs (Glynn, 1996). Blank (1987) recon¬ 
structed the nucleus of Symbiodinium kawagutii and found 
that the chromo.somes of this dinoflagcllate could not be 
paired either by size, appearance, or distribution. This cy¬ 
tological result led to the speculation that the vegetative 
(coccoid) cells of Symbiodinium may be haploid (Blank, 
1987). To date, the molecular genetic data necessary to 
establish the ploidy of these symbiotic dinollagcllates, or 
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any dinollagellate, are lacking. We resolve this lingering 
question by using inicrosatellites to elucidate the ploidy of 
Symbiodiniwu. 

Microsatellites are simple, tandemly repeated DNA se¬ 
quences (reviewed in Chambers and MacAvoy, 2000) that 
are distributed abundantly in the genomes of virtually all 
organisms (Bennett, 2000). These single-locus, multiallelic, 
codominant segments of DNA are highly versatile and 
accessible markers for population genetic studies. However, 
micro.satellites have also been successfully applied as tools 
in determining an organism’s ploidy. For example, recovery 
of two alleles from a single locus in a single individual or 
isolate demonstrates a diploid nuclear phase. On the other 
hand, if a single allele is recovered from each of multiple 
loci, and the result is replicated over multiple individuals or 
isolates, the data suggest a haploid nuclear phase. This 
rationale has been applied to numerous organisms, includ¬ 
ing a parasitic protozoan, Tiypanosoma cntzi (Oliveira et 
al., 1998); two parasitic fungi. Maguapoithe grisea (Bron- 
dani et ai, 2000) ^ind Aspergillus finni gat us {Bavi-DcVdbQssQ 
et ai, 1998); a diatom, Ditylum briglitwellii (Rynearson and 
Armbrust, 2000); a bryophyte, Polytricluuu foniiosuni (van 
der Velde et ai, 2001); and a false spider mite, Brevipalpiis 
pboenicis (Weeks et ai, 2001). In this investigation, we 
apply the same strategy to members of the genus Symbio- 
diniwjh Populations of Synibiodinium harbored by the gor- 
gonians Plexaiira kuna and Pseudopterogorgia elisabethae, 
as well as a series of Symbiodiiiiuin cultures, were screened 
with five polymoiphic microsatellites. The results confirm 
that members of Syuibiodinium are haploid in the vegetative 
life stage. 

Materials and Methods 

Biological materials and nucleic acid isolation 

Samples of Plexaura kuna were collected from colonies 
at depths of 1-7 m at 10 sites in the San Bias Islands, 
Republic of Panama {n = 142); two sites in the Florida 
Keys {n = 12); and three sites in the Bahamas (n = 6). For 
Pseudopterogorgia elisabetliae, samples were collected 
from 575 colonies at depths of 10-27 m at 12 sites in the 
Bahamas (n = 43-50 per site). The 12 P. elisabetluie sites 
were Sweetings Cay, Gorda Rock, Little San Salvador, Cat 
Island, Rum Cay, Hog Cay, and two sites each at Abacos, 
Eleuthera, and San Salvador Islands. Immediately after col¬ 
lection. samples were preserved in either 95% ethanol or 
salt-saturated DMSO (Seutin et ai, 1991) or were frozen in 
a liquid nitrogen vapor shipper. Total nucleic acids were 
extracted and quantified from branch pieces (about 2 cm) of 
P. kuna according to the methods of Coffroth et ai (1992). 
Nucleic acids were extracted from branch pieces (about 3 
cm) of P. elisabetluie by first grinding the tissue in STE 
buffer (0.05 M Tris-HCl (pH = 8.0), ().l M EDTA, 0.1 M 
NaCl, 0.2% SDS) followed by a modified extraction proto¬ 


col using the Prep-A-Gene DNA Extractioii Kit (Bio-Rad 
Laboratories, Hercules, CA) (Shearer and Coffroth, un- 
pubL). These extraction methods extracted nucleic acids 
from both the gorgonian hosts and from their Synibiodinium 
populations in hospite. 

Synibiodinium cultures, isolated from a range of inverte¬ 
brate hosts and geographic locations (Table 1), were main¬ 
tained as described in Santos et ai (2001). Total nucleic 
acids were extracted and quantified from fresh algal cultures 
(approximately 5 X 10^ cells) as described above for the P. 
kuna samples. After extraction, all nucleic acid samples 
were stored at —20 °C. 

Phylogenetic identification o/Synibiodinium populations 
and cultures 

The phylogenetic identity of the Synibiodinium popula¬ 
tions harbored by P. kuna and P. elisabetluie, as well as the 
Synibiodinium cultures, was determined by length hetero- 
plasmy in domain V of chloroplast large subunit (cp23S) 
rDNA as described in Santos et ai (2003). This technique 
rapidly genotypes Synibiodinium isolates and places them 
into a phylogenetic framework based on cp23S-rDNA (San¬ 
tos et ai, 2003). 

Construction o/Synibiodinium microsatellite libraries 

Nucleic acids used in the construction of the microsatel- 
lite libraries for P. kuna and P. elisabetluie Synibiodinium 
populations were obtained by different procedures. Symbio- 
diniiim was isolated from samples of P. kuna collected from 
around the Caribbean and purified of host tissue (described 
in Santos et ai, 2001) before the nucleic acids were ex¬ 
tracted and quantified as described above for P. kuna sam¬ 
ples. On the other hand, nucleic acids for the Synibiodinium 
microsatellite library from P. elisabetluie were extracted 
from colonies collected at Sweetings Cay in = 2) and San 
Salvador (n — 4) as described above. These preparations 
contained nucleic acids from both P. elisabetluie and their 
Synibiodinium populations in hospite, with no effort being 
made to enrich for Symbiodiniwu nucleic acids. After ex¬ 
traction, nucleic acids were pooled according to host source 
and used to construct the microsatellite libraries. The librar¬ 
ies were constmcted and screened for microsatellites with 
dinucleotide repeats as described in Cioti and Bruford 
(1998). 

Screening of microsatellite loci polymerase chain reaction 
primers 

To ensure that the microsatellite loci were not part of the 
host gorgonian genome, PCR amplifications were carried 
out on nucleic acids extracted from planulae of P. kuna and 
P. elisabetluie, which are asymbiotic upon release from the 
maternal colony (Kinzie, 1974; Coffroth et ai, 2001; see 
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Table I 


Information about the Synibiodiniiim cultures that amplifieJ with at 
least one of the two microsatelfites primer sets designetf for the 
Symbiodinium populations Pseudopterogorgia elisabethae 


Host organism 

Culture name 

Collection location 

Aiptasia pallUla 

FLAp#2 

Florida Keys 

A. pallUIa 

Fl.Ap#2 lOAB 

Florida Keys 

A. palliiJa 

FLAp#3 lOAB 

Florida Keys 

A. pallida 

FLAp#4 lOAB 

Florida Keys 

Aiptasia pulchella 

OkAp#l 

Okinawa, Japan 

A. pulchella 

(4kAp#9 

Okinawa, Japan 

A- pulchella 

OkAp#l() 

Okinawa, Japan 

Briareum ashestinum 



(polyp)'' 

#498^* 

Florida Keys 

B. ashestinum (polyp 

#579’' 

Florida Keys 

B. ashestinum (polypP 

#595^^ 

Florida Keys 

B. ashestinum (polyp)** 

#1135” 

Florida Keys 

B. ashestinum (polyp)** 

#1140^ 

Florida Keys 

B. ashestinum (polyp)** 

#1246^ 

Florida Keys 

B. ashestinum (polyp)** 

#1319” 

Florida Keys 

B. ashestinum (polyp)** 

#1385*^ 

Florida Keys 

B. ashestinum (polyp)** 

#1394^ 

Florida Keys 

B. ashestinum (polyp)** 

#1509'^ 

Fhirida Keys 

B. ashestinum (polyp)** 

#1510'^ 

Florida Keys 

B. ashestinum (polyp)** 

#1(^69'^ 

Florida Keys 

B. ashestinum (polyp)** 

#1735'^ 

Florida Keys 

B. ashestinum (polyp)** 

#1902^ 

Florida Keys 

B. ashestinum (polyp }"* 

#2053^ 

Florida Keys 

B. ashestinum (polyp)** 

#2080^' 

Florida Keys 

Plexaura kuna 

Pkl3 

Florida Keys 

P. kuna 

Pkl4 

Florida Keys 

P. kuna (polyp)** 

Pk2()5 

San Bias Islands, Panama 

P. kuna (polyp)** 

Pk2()6^ 

San Bias Islands. Panama 

P. kuna (polyp)** 

Pk208^ 

San Bias Islands. Panama 

P. kuna (polyp)** 

Pk215 

San Bias Islands, Panama 

P. kuna (polyp)** 

Pk216 

San Bias Islands, Panama 

P. kuna (polyp)^ 

Pk225'* 

San Bias Islands, Panama 

P. kuna (polyp)** 

Pk226^ 

San Bias Islands. Panama 

P. kuna 

Pk8()l 

Florida Keys 

P. kuna 

Pk807 

Florida Keys 

P. kuna (polyp)** 

Pk7()2‘^ 

San Bias Islands, Panama 

P. kuna (polyp)** 

Pk703^ 

San Bias Islands. Panama 

P. kuna (polyp)^ 

Pk7n4^ 

San Bias Islands, Panama 

P. kuna (polyp)** 

Pk705'* 

San Bias Islands, Panama 

P. kuna (polyp )"* 

Pk7()6'" 

San Bias Islands, Panama 

P. kuna (polyp)** 

Pk7n7^’ 

San Bias Islands, Panama 

PtfciUoptna damieornis 

Pd 

Hawaii 

Porites evennanni 

Pc 

Hawaii 

Pseiidt >pter( ^i'o rg ia 



elisabethae 

SSPc 

Bahamas 

Unknown anemone 

Cane 

Hawaii 

Zoanthus paeificus 

Zp 

ttawaii 


^ A polyp is defined as a newly settled and nietamcirphosed planula. 
^Cultures that were started from a single dinoflagellate cell. 


below). Nucleic acid.s were extracted froiii the plaiuilae as 
described above for samples of P. kuna. To lest for the 
presence of template nucleic acids {i.e., cnidarian host 
DNA), FCR amplifications were conducted using the uni¬ 


versal nI8S-rDNA primer set ss5 and ss3 (Rowan and 
Powers, 1991). The asynibiotic nature of the planulae was 
then confirmed by the absence of a PCR product when the 
zooxanthella-biascd nl8S-rDNA primers 5s5 and ss3z were 
used (Rowan and Powers, 1991). Both sets of reactions 
were carried out under the conditions described by Rowan 
and Powers (1991). Details of the five Symbiodinium mi- 
crosalellite loci, including the sequences of the PCR primers 
and GenBank accession numbers, are given in Table 2. 

Micwsatellite amplifications and allele detection for 
Symbiodinium populations 

PCR reactions for each Symbiodinium microsatellite lo¬ 
cus were performed in 10-/xl volumes containing 10 mA/ 
Tris-HCl (pH 8.3), 50 niA/ KCl, 0.001% gelatin, 200 /xA/ 
dNTPs, 1 U Tacj polymerase, approximately 10 ng of tem¬ 
plate DNA, and MgCH and primers at concentrations de¬ 
tailed in Table 2. Thermocycling conditions were as fol¬ 
lows: initial denaturing step of 94°C for 2-3 min, 35-39 
cycles of 94°C for 30-45 s, annealing temperature of X° 
(see Table 2) for 30-45 s, 72®C for 30 s, and a final 
extension of 72°C for 3-5 min. 

Most Symbiodinium micro.satellite alleles from P. kuna 
were separated in 2% 0.5X Tris-borate (TBE) agarose gels 
and visualized by ethidium bromide staining. In addition, 
some Symbiodinium microsatellite alleles from P. kuna 
were amplified in the presence of ^'^P-ATP (Sambrook et ai. 
1989), separated in 6% polyacrylamide / 0.5X TBE gels 
under denaturing (7 A/ urea) conditions, and visualized by 
exposure to X-ray film. The sizes of Symbiodinium micro¬ 
satellite alleles from P. kuna were determined either with 
a 100-bp DNA ladder (MBl Fermentas, Hanover, MD) 
for 2% agarose gels, or for polyacrylamide gels, with a 
^'^P-ATP puC19 plasmid sequencing reaction serving as a 
nucleotide size-ladder. 

For P. elisabethae, Symbiodinium microsatellite alleles 
were labeled by incoiporalion of a 5'-lRD800 M13 Forward 
primer (see Table 2), and were separated in 25-cm-long, 
0.25-mm-thick, 6.5% Long Ranger (FMC Bioproducts, 
Rockland, ME)/0.5X TBE gels under denaturing (7 M urea) 
conditions. Gel electrophoresis was performed at 1500 V, 
40 W, 40 mA, 50 ®C, and the default scan speed, with 
Ll-COR's NEN Global 1R2 DNA sequencer system 
(Ll-COR Biotechnology Division, Lincoln, NE). The sizes 
of Symbiodinium microsatellite alleles from P. elisabethae 
were determined with the fragment analysis program Gene 
ImaglR v3.55 (Scanalytics Inc, Fairfax, VA) using DNA 
ladders (ladder 97-147 bp or 172-272 bp, rungs at 25-bp 
increments, ladders loaded every fifth lane) as size refer¬ 
ences. 

In addition, Symbiodinium microsatellitc-loci primer-sets 
for each host species were tested on the Symbiodinium 
populations of the other host species. For the host species 
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I'nble 2 


Characterization and PCR amplification conditions of the five Symbiodinium clade B microsatellite loci: size refers to the predicted PCR product size 
w ith original seciiienced done as template 


Microsaieltite 

locus name 

Motif type 

Sequence of 
microsatelliie motif 

PCR primer sequences 

Size 

(bp) 

cAr.?'*’ 

GT, mlemipled pure 

(GTi.GATtGTCGC 

(GT),3GAT(GTL 

CAl '.7UP: 5'-GCAGTGTGn'ACCrTAAAGTGCGG-3' 

CA 1 '.7DN: 5'-GTTTAGTTGAGCAACTTCCGAG-3' 

345 

GA2.8^'’ 

CA, interrupted pure 

(CAt^CGtCAL, 

GA2.81JP; 5'-TCGACTCTTGGCGCAAAATG-3' 

GA2.8DN: 5'-GGGATGAAACACTGGGATAATCCAG-3' 

174 

GA4.84'‘-'’ 

CA, interrupted complex 

(CAGAmCAlCGA 

A(CAL(CG).CAGA 

(CA),3CG(CA)3 

(GA(CA).),GACG 

(CA), 

GA4.S4UP: 5'-GATCAAACTCTTGTGATGAC-3' 

GA4.84DN; 5'-GTCAGATTGTCA4CAAAGACrGC-3’ 

175 

CA4.86'''^ 

CA, interrupted complex 

(CACC)CA 

(CACOaCA)^ 

CA4.86R; 5'-GCCITCAATGCAATCACCTT-3' 

C A4.86L': 5 '-GG A AT IGGCCATCCCTCT AT-3' 

193 

CA6.38'--‘^ 

CA, interrupted pure 

(CA).,cg((:a)4 

C A6.38R: 5 '-CA A AG A AT ATTCGGGGGTC A-3 ’ 

CA6.38I3: 5'-AGTTGATACGCCGGA rGTGT-3' 

112 


^ Microsaielliie loci isolated from Plexaura kuna Symbiodinium populations in hospite; sequences deposited in GenBank under accession numbers 
AF549186-AF549I88. 

^ 1.5 xnM MgCE. 3 pniol of each primer, 60 annealing temperature. 

^ Microsalellite loci isolated from Pseudopterogorgia elisabethae Symbiodinium populations in hospite; sequences deposited in GenBank under accession 
numbers AF474165 and AF474169. 

'*2.5 mA/ MgCln. 0.2 pmol CA4.86R, 0.18 pmol CA4.86LM13, 0.02 pmol 5'-lRD800 M13 Forward primer, 50 annealing temperature. 

1.5 xxh\f MgCU, 3 pmol CA6.38R, 2 pmol CA6.38FM13. 0.2 pmol 5'-IRD800 M13 Forward primer, 56 X annealing temperature. 

* Primers with M13 Forward sequence (5'-CACGACGTTGTAAAACGAC-3') added for visualization on Ll-COR's NEN® Global 1R2 DNA Sequencer 
System (LI-COR Biotechnology Division, Lincoln, NE). See text for details. 


comparisons, PCR reaction conditions and detection of 
Symbiodinium microsatellite alleles were as previously de¬ 
scribed, 

A total of 45 Symbiodinium cultures were also .screened 
with the five Symbiodinium microsatellite loci primer sets 
using the methods described above. 

Statistical testing of microsatellite data 

Exact probabilities (analogous to Fisher's exact lest for 
2X2 contingency tables) for ihe observed genotype fre¬ 
quencies and for those expected under Hardy-Weinberg 
equilibrium were used to determine whether populations of 
Symbiodinium harbored by P. elisabethae in the Bahamas 
deviate from Hardy-Weinberg equilibrium. Significance 
values were calculated with the computer program BIO- 
SYS-2 (Sw'offord and Selander. 1981; modified by William 
C. Black. Deparlment of Microbiology, Colorado State Uni¬ 
versity; available at: ftp://lamar.colostate.edu/pubAvcb4/). 

Results 

Phylogenetic identification qf Symbiodinium populations 
and cultures 

The numerically dominant Symbiodinium in association 
with Plexaura kuna and Pseudopterogorgia elisabethae be¬ 


long to Symbiodinium clade B {sensu Rowan and Powers, 
1991), and more specitically, to Symbiodinium B184 (sensu 
Santos et ah, 2003). The Symbiodinium ciillures also belong 
to clade B, with represenlatives from Symbiodinium B184, 
B211, and B223 (sensu Santos et al., 2003; Table 3). 

Genomic location of the microsatellite loci 

Nucleic acids extracted from planulae of P. kuna and P. 
elisabethae failed to produce an amp!icon with ihe zooxan- 
thella-biased nl8S-rDNA PCR primers ss5 and ss3z. On the 
other hand, an arnplicon of approximalely 1800 bp was 
generated from the same templates using the universal 
nl8S-rDNA PCR primers ss5 and ss3 (dala not shown). 
These results confirmed that the planulae were asymhiolic. 
Nucleic acids extracted from planulae and then used as 
templates for the microsatellite loci primer sets failed to 
produce PCR amplicons, whereas nucleic acids from adult 
tissue that contained Symbiodinium or nucleic acids isolated 
from Symbiodinium cultures produced PCR amplicons (.see 
below) of the correct size (given in Table 2). Taken to¬ 
gether. the data demonstrate that the microsatellite loci are 
located within the genome o\ Symbiodinium, and not that of 
the gorgonian hosts. 
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Table 3 


Micmsateiiite ami chtoroplast genotypes for the Symbiodinium cultures: a microsateltite 
allele sizes at loci CATcS‘6 ami CA6.3S 

genotype is defined as a 

uniipte comhination of microsatellite 

SYmhioJinium microsalellite allele 




size (bp) 

Sxiuhioiliniuin chloroplast 




large subunit rDNA 



Locus CA4.86 Locus CA6.38 

(cp23S-rDNA) genotype'" 

Culture name 

Host organism 

179 100 

B184 

Cane 

Unknown anemone 


B184 

OkAp#l 

Aiptasia pidchella 


B184 

OkAp#l() 

A. pidchella 


B184 

OkAp#9 

A. pidchella 


B184 

Pd 

Pocillopora daiuiconiis 


B184 

Pe 

Foiites evernianni 


B184 

Zp 

Zoantluis pacificus 

179 102 

B184 

FLAp#2 

Aiptasia pallida 


B184 

FLAp#2 lOAB 

A. pallida 


B184 

FLAp#3 lOAB 

A. pallida 


B184 

FLAp#4 !()AB 

A. pallida 

183 100 

B184 

Pk205 

Plexaura kinia (polyp)" 


B184 

Pk208^ 

P. kuna (polyp)"" 


B184 

Pk2LS 

P. kuna (polyp)"" 


B184 

Pk216 

P. kuna (polyp)"" 


B184 

Pk226‘^ 

P. kuna (polyp)"" 

187 ‘nulF’ 

B223 

#498^ 

Briareiun asbestinum (polyp)" 


B223 

#579'’ 

B. ashestiniuii (polyp)" 


B223 

#1385'’ 

B. asbestinum (polyp)" 

191 104 

B184 

#595*^ 

B. asbestinum (polyp)" 


B184 

#1246'’ 

B. asbestinum (polyp)" 


B184 

#1394'’ 

B. asbestinum (polyp)" 

191 106 

B184 

#1 135'’ 

B. asbestinum (polyp)" 


B184 

#1140'’ 

B. asbestinum (polyp)" 


B184 

#1319'’ 

B. asbestinum (polyp)" 


B184 

#1509'’ 

B. asbestinum (polyp)" 


B184 

#1510'’ 

B. asbestinum (polyp)" 


B184 

#1669'’ 

B. asbestinum (polyp)" 


B1S4 

#1735'’ 

B. asbestinum (polyp)" 


B184 

#1902'’ 

B. asbestinum (polyp)" 


B184 

#2053'’ 

B. asbestinum (polyp)" 


B184 

#2080'’ 

B. asbestinum (polyp)"' 


B1S4 

Pk225'’ 

P. kuna (polyp)"' 


B184 

Pk706'’ 

P. kuna (polyp)" 


B184 

Pk7n7'’ 

P. kuna (polyp)" 

193 98 

B2I 1 

Pk702'’ 

P. kuna (polyp)" 


B21 1 

Pk703'‘ 

P. kuna (polyp)" 

193 100 

B184 

Pkl3 

P. kuna 


B184 

Pkl4 

P. kuna 


B184 

PkSOl 

P. kuna 


B184 

Pk807 

P. kuna 

193 104 

B184 

Pk206'’ 

P. kuna (polyp)" 


B184 

Pk704'’ 

P. kuna (polyp)" 


BI84 

Pk705'’ 

P. kuna (polvp)" 

193 112 

B184 

SSPe 

Pseudopterogorgia elisabetliae 


*' Scnsu Samos et al. (2003), 

Cultures iliai were started from a single dinotlageilaie eell. 

A polyp is defined as a newly settled and nieiamorphosed planuia. 


Microsatelliles Symbiodiiiiuin within a host species 

For I he Synihiodinimn populations of l\ kuna, four 
unicjuc alleles were identilied at locus CAT.?, three al 


GA2.S, and three at GA4.84. The sizes of the alleles at each 
locus were CAT.7 (approximately 360, 380, 400, 420 bp), 
GA2.8 (approximately 180. 190, 2()() bp), and GA4.84 (ap¬ 
proximately 280, 300, 310 bp). In 147 cases, only a single 
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allele per locus could be delected from a P. kuna colony, but 
the Syinbiodinium populations from 13 colonies produced 
two distinct alleles at one or more loci. Thus, 8.1% of the 
Symhiodinium populations sampled from P. kuna individu¬ 
als produced more than a single allele at any microsatellite 
locus. 

When the Synihiadiniuni populations of P. elisahethae 
were screened at loci CA4.86 and CA6.38, 8 and 10 unique 
alleles, respectively, were identified from 568 of the 575 
colonies that were examined (nucleic acids from seven 
colonies failed to amplify at either locus and were excluded 
from further analysis). Allele sizes for loci CA4.86 and 
CA6.38 ranged between 185-207 bp and 96-122 bp, re¬ 
spectively (Fig. 1). In most of these cases, only a single 
allele per locus could be detected from the Symbiodinium 
population of a P, eUsabetbae colony (Fig. 2). In 25 cases 
(4.4%), the Symbiodinium population from a P, eUsabetbae 
colony possessed two distinct alleles at one or more loci 
(example in Fig. 2). This pattern is similar to that observed 
for Symbiodinium microsatellites isolated and amplified 
from P. kuna. 


Symbiodinium microsatellites between host sjyecies 

Amplifications using the Symbiodinium microsatellite 
loci primer sets for one host species produced mixed results 
when applied to the Symbiodinium populations of the other 
host species. In many cases, the primers failed to produce an 
amplicon. However, when PCR amplification did occur, 
only a single allele was detected per locus (see hg. 4 of 
Santos et ak, 2001, for examples). 

Microsatellites from Symbiodinium cultures 

The Symbiodinium microsatellite loci CAF.7, GA2.8, 
and GA4.84 were not detected in any of the Symbiodinium 
cultures. However, the 45 cultures did amplify with at least 
one of the two Symbiodinium microsatellites primer sets 
isolated from P. eUsabetbae (Table 2). At loci CA4.86 and 
CA6.38, 5 and 6 alleles, respectively, were identified (Table 
3, Fig. 1). Allele sizes for loci CA4.86 and CA6.38 ranged 
between 179-193 bp and 98-1 12 bp, respectively (Fig. 1). 
In each case, however, only a single allele was detected per 
locus. 
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Figure 1. Size disiribuhon of alleles al two microsaielliie loci; (a) 
CA4.86 and (b) CA6.38. Solid bars. SynihiocHtiinm populations of Pseu- 
doptero^orgia elisahethae {in hospite); diagonal sirip bars, Synihiadiniuni 
cultures. Values for P. elisahethae are derived from all samples, include 
those that possessed two distinct alleles or a null allele at one or more loci. 
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Test for deviations from Hardy-Weinberg equilibrium 

The observed genotype frequencies of the Symbiodinium 
populations inhabiting Pseudopterogorgia eUsabetbae from 
most of the 12 sites in the Bahamas were significantly 
different from those expected under Hardy-Weinberg equi¬ 
librium (Table 4). An excess of "homozygote" (genotypes 
with a single allele per locus) genotypes w'as present in the 
populations compared to the expected number based on 
allelic frequencies. In four cases, a single monomoi*phic 
allele w^as recovered from the population (Table 4). 

Discussion 

Microsatellite loci isolated IVom members of Symbio¬ 
dinium clade B were amplified by PCR to assess ploidy in 
these symbiotic dinofiagellates. From most Symbiodinium 
populations harbored by Ple.xaura kuna and Pseudoptero¬ 
gorgia eUsabetbae (690 out of 728 samples: 95% of the 
cases), only a single allele was recovered per locus, and all 
45 Symbiodinium cultures possessed a single allele per 
locus. Furthermore, an excess of "homozygous," or single 
allele per locus, genotypes in the Symbiodinium populations 
of P. eUsabetbae violated Hardy-Weinberg equilibrium. 
Nonconformity to the prediction of Hardy-Weinberg equi¬ 
librium indicates that one or more of its a.ssumptions are not 
met in the population. The assumptions include (1) the 
organism is diploid; (2) reproduction is sexual; (3) the 
population is infinitely large: (4) mating occurs randomly; 
(5) generations are nonoveiiapping; and (6) the population 
is free of genetic drift, migration, mutation, and natural 
selection (review^ed in Hartl and Clark. 1989). The Symbio- 
ditiium populations of P. eUsabetbae may violate many of 
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Figure 2. Symhioditiuim populalions from Fseudopterogorgki elisid->etlhie: Polyacrylamide gel electro¬ 
phoresis analysis of microsatellite alleles from two loci: (a) CA4.86 and (h) CA6.38. Number above lane 
represents individual samples of P. elisabcthae. L represents DNA ladder lanes. Vertical arrows denote samples 
in which more than a single allele was amplified per locus. 


these assumptions, but the predominance of "homozygous" 
genotypes in tlie Syuibiodiniuni populations harbored by P. 
kuna, and in all Syuihiodiuuuu cultures, strongly suggests 
that at least the first assumption is being violated. If Syui- 
hiodiniuiu were diploid in tlie vegetative life stage, more 
"heterozygous" genotypes should have been sampled. A 
possible exphtnalion for the lack of heterozygous genotypes 
is that the Symhiodiuium microsatellite loci are located in an 
organellai (/.r\. mitocliondrial or chloroplast) genome. If so, 
a similar pattern of a single allele per locus would be 
observed. Although a mierosatellite locus has been de¬ 
scribed from the chloroplast genome of a free-living 


dinoflagellate, lieterocapsa tric/uetra (Zhang ef ok, 1999), it 
is highly unlikely that all five Syuihiodiiuinu mierosatellite 
loci were isolated from the organellar genomes. Taken 
together, our results demonstrate that members of Symbio- 
diuiwu clade B, both cultured and in Iwspite, possess a 
haploid nuclear phase. Moreover, since the genus is mono- 
phyletic (see, for e.xample. Rowan and Powers, 1992; La- 
Jeunesse, 2001; Santos ct ai, 2002), our findings coiTobo- 
rate Blank's (1987) speculation that haploidy e.xisls in S. 
kawaputii and extend it to all members of Syiubiodiuiwu. 

If the vegetative life stage of Syuibiodiuiiini is haploid, 
why was more than a single allele sometimes recovered per 
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Table 4 


Significant deviation of ohsen'ed genotype freqnencies from those 
expected under Hardy-Weinberg e(]nilihrinnj for Symbiodinium 
populations inhabiting P^.eudopterogorgia clisabethae in the Bahamas 


Naine/localion of Psendopterogorgia 
clisabethae populations (sample size) 

P values 

Locus CA4.86 

Locus CA6.38 

Sweetings Cay (44) 

<0.001 

<0.001 

Gorda Rock (47) 

<0.001 

<0.001 

Abacos Shallow (50) 

<0.001 

<0.001 

Abacos Deep (48) 

<0.001 

<0.(J01 

South Hampton Reef (47) 

0.011 

0.032 

Little San Salvador (48) 

<0.001 

<0.001 

East End Point of Eleuthera (44) 

<0.()()1 

<0.()0t 

Cat Island (45) 

<0.001 

MA 

San Salvador, Riding Rock (43) 

MA 

MA 

San Salvador. Pillar Reef (50) 

<0.001 

MA 

Rum Cay (46) 

<0.(J01 

<0.001 

Hog Cay (46) 

<().()01 

<0.001 


Values were calculated using exact probabilities (analogous to Fisher's 
exact test tor 2 X 2 contingency tables) using the computer program 
BlOSYS-2 (Swofford and Selander. 1981; modified by William C. Black, 
Department of Microbiology, Colorado State University; available at: 
ftp://lamar.colostate.edu/pub/wcb4/). Samples that possessed a “null.” or 
absent, allele at a locus were excluded from the analysis. MA designates a 
monomorphic allele recovered from the population. 


locus in hospite? Among the samples of P, kuna and P, 
elisahethae, about 8.1 % and respectively, gave this 

result. In these cases, the recovery of more than a single 
allele at a microsatcllite locus suggests that the colonies 
harbored at least two genotypes of symbiotic algae. It has 
been demonstrated, using a combination of culturing and 
molecular techniques, that P. kuna and P. elisahethae can 
harbor more than a single genotype of SynihiocHnium simul¬ 
taneously (Goulet and Coffroth, 1997, 2003: Coffroth et ak, 
2001; Santos et ak, 2001). Our data support this conclusion 
while suggesting lhat it is uncommon for colonies of these 
gorgonians to harbor more than one Syinbiodiniwn geno¬ 
type, at least at detectable concentrations. For these uncom¬ 
mon colonies, it remains to be determined whether the 
additional genotypes represent relicts of the initial symbiont 
uptake by the newly settled asymbiotic planulae or reflect 
populations that are secondarily acquired and contribute to 
the host in other ways. 

Phylogenetic support for haplouly in Symbiodinium and 
other dinoflageJlates 

Our conclusion that Sytnbiodiniuni is haploid is consistent 
with data from most of the other dinollagellates. The life 
cycle of nearly all dinollagellates examined to date is dom¬ 
inated by asexual reproduction of haploid vegetative cells 
(Pficster and Anderson, 1987). This observation, coupled 


with the monophyly of dinollagellates (Saldarriaga et ak, 
2001), suggests that all members of the dinollagellates 
should possess a vegetative haploid nuclear phase. Interest¬ 
ingly, the vegetative cells of the red-tide dinollagellale 
genus Noctiluca are thought to be diploid (Zingmark, 1970; 
Pfiester and Anderson, 1987). In Noctiluca, the first divi¬ 
sions of the gamete mother-cell nucleus are believed to be 
meiotic (Zingmark, 1970), which would imply dipioidy in 
these dinoflagellates. Recently, the conclusion of dipioidy in 
Noctiluca has been challenged (Schnepf and Drebes, 1993). 
but no definilive data have been presented to establish the 
ploidy of these dinollagellates. We hypothesize, based on 
our knowledge of other dinollagellates. that Noctiluca spp. 
possess a vegetative haploid nuclear phase. Analyses of 
microsatellile loci, such as we have done here for Sytnbio- 
diniiun, would be one way to test this hypothesis and settle 
the question of ploidy in Noctiluca. 

A haploid nuclear phase in the dinollagellates is consis¬ 
tent with that of their closest relatives. The Apicomplexa, 
obligate intracellular parasites of many vertebrate and in¬ 
vertebrate hosts, are thought to have evolved from, or 
shared a common ancestor with, the dinollagellates 
(Woliers, 1991: Cavalier-Smith, 1993) aboul 395-929 Mya 
(Escalante and Ayala, 1995). The apicomplexan Plasmo¬ 
dium falciparum, one of the causative agents of human 
malaria, is haploid in its human host and only briefly diploid 
in its mosquito vector (Campbell, 1993; Conway et ak, 
1999). Other apicomplexans, such as Ctyptosporidium par- 
vum and Toxoplasma gondii, also possess a haploid nuclear 
phase (Costa et ak, 1997; Feng et ak, 2002). Given the close 
evolutionary relationship between the two groups, the an¬ 
cestral state in the progenitor of the apicomplexans and 
dinollagellates was probably haploidy. 

Evidence for fine-scale specificity in associations betweeti 
host and Symbiodinium 

Surprisingly, primer sets designed for amplification of the 
Symbiodinium microsatellile loci in one host species pro¬ 
duced mixed results when applied to the Symbiodinium 
populations of the other host species or to Symbiodinium 
cultures. For example, primers designed for the symbiont 
populations of P. kuna were not very successful in ampli¬ 
fying the Symbiodinium populations of P. elisabethae col¬ 
onies or algal cultures derived from a variety of hosts 
(Santos et ak, 2001; unpubl. data). Typically, the utility of 
a microsatellile system {i.e., microsalellite primer sets) de¬ 
creases as the phylogenetic distance between the samples 
being screened increases (Schlotterer, 1998). In these ex¬ 
periments, however, most (768 out of 773; 99.4%) of the 
samples belonged to a single lineage, Symbiodinium B184. 
Therefore, the microsatellile primer sets should have 
worked on all members of the group. This failure to amplify 
alleles from closely related, but non-focal. Symbiodinium 
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samples is probably due to mutational changes in the flank¬ 
ing regions of the microsatellite airay. Mutations in these 
regions can lead to primer-template mismatch and thus to 
inhibition of the PCR reaction. In support of this hypothesis, 
we have sequenced alleles from loci CA4.86 and CA6.38 
and found mutations, such as nucleotide substitutions and 
insertion-deletions (indels), in the microsatellite flanking 
regions from members of SyDihiodifiiwii B184, B211, and 
B223 (the evolution of Synihiodiniu)}! microsatellites will be 
discussed in a subsequent paper). 

Although microsatellite alleles were not always recov¬ 
ered in the host species comparisons, an important conclu¬ 
sion can be drawn from these data. The specificity exhibited 
by the different microsatellite primer sets to the population 
of Sy))d}i()difiiit})i from which they were designed and the 
consistent presence of “niilf" or absent, alleles in the other 
host species suggest that the two Syiuhiodiiiiinn B184 pop¬ 
ulations are genetically distinct from each another and spe¬ 
cific to a given host species. The genetic differences are 
probably spread across the Sy)uhiodiniuf}i genome, but at 
the minimum they are confined to mutations in the flanking 
regions of the microsatellite loci. Unfortunately, internal 
transcribed spacer (ITS) sequences—one of the most useful 
genetic markers for identifying Syiiddodiiuiifu types (LaJeu- 
nesse, 2001)—are identical, or nearly so, in the two popu¬ 
lations (Santos et cd., 2001): thus other genetic markers are 
needed to elucidate the relationship between them. Never¬ 
theless, P, kuna and P. elisahethae appear to associate 
preferentially with genetically distinct Syudyiodiiiiuui B184 
populations, which provides e\ idence for fine-scale hos\- 
Syudyiodiuiiini specificity in these gorgonian species. 

MicrosatcUiies and Symbiodinium diversity 

Consistent with other studies (Schoenberg and Trench, 
1980; Colley, 1984; Goulet and Coffroth, 1997, 2003; Bail- 
lie et ai, 1998, 2000; Belda-Baillie et a!., 1999), our mic- 
rosalellite data suggest an enormous amount of genotypic 
diversity within Synihiodiniuiu, as illustrated by the follow¬ 
ing example. At loci CA4.86 and CA6.38, a total of 10 and 
12 unique alleles, respectively, were recovered from sam¬ 
ples belonging to Synihiodiniiiin B184. Pairing alleles from 
each locus, under the assumption that there are no restric¬ 
tions against particular combinations of alleles, generates 
120 unique genotypes of Synddodiiiiinu B184. However, we 
feel that this is a conservative estimate for several reasons. 
First, a minimal number of microsatellite loci are being 
employed. Data from other polymorphic microsatellite loci 
would distinguish more genotypes within Symbiodiuiiun 
B184. Second, some alleles for CA4.86 and CA6.38 are 
missing from the data set because they have not yet been 
.sampled (Fig. 1). The inclusion of any of these alleles would 
generate up to 210 unique genotypes of Symbiodinium 
B184. Third, the Symbiodinium BI84 populations of hosts 


such as P. kuna possess “null" alleles at these loci, suggest¬ 
ing an additional level of diversity within the group (see 
above). Last it is extremely unlikely that this high level of 
genotypic diversity is confined to Symbiodinium B184. 
Thus, microsatellites will doubtlessly uncover high levels of 
genotypic diversity in most, if not all, of the 16 Symbio¬ 
dinium lineages recognized in cp23S-rDNA phylogenies 
(Santos et ai, 2002, 2003). 

Evidence for recofidnnation in Symbiodinium 

The finding that vegetative cells of Syiiibiodinium possess 
a haploid nuclear phase does not preclude recombination 
within the life cycle of these symbiotic dinoflagellates. For 
example, other haploid organisms maintain some form of 
recombination during their life cycle, including the green 
alga Chlamydonionas and members of the Acrasiomycota 
(cellular slime molds), the Bryophyta (mosses), the Ptero- 
phyta (ferns), the Apicomplexa, and the Dinophyceae 
(dinoflagellates) (Pfiester and Anderson, 1987; Campbell, 
1993). In fad, the high allelic variability observed for 
allozymes (Schoenberg and Trench, 1980; Baillie et a!., 
1998; Belda-Baillie et aL, 1999), random-amplified-poly- 
morphic DNA (RAPDs) (Belda-Baillie et ai, 1999; Baillie 
et ai, 2000), and DNA fingeiprints (Goulet and Coffroth, 
1997, 2003) suggests extensive recombination in Symbio¬ 
dinium (Baillie et ai, 2000; review-ed in LaJeunesse, 2001). 
This evidence for recombination, taken together with our 
finding of haploidy, lends strong support to Symbiodinium 
life cycle (a), as proposed by Fitt and Trench (1983). How¬ 
ever, questions pertaining to recombination in these enig¬ 
matic dinoflagellates, such as the factors that induce it and 
whether it occurs inside or outside a host, remain to be 
answ'ered. 
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